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The ground state of a quantum spin chain is a natural playground for investigating cor- 
relations. Nevertheless, not all correlations are genuinely of quantum nature. Here we 
review the recent progress to quantify the 'quantumness' of the correlations throughout 
the phase diagram of quantum spin systems. Focusing to one spatial dimension, we dis- 
cuss the behavior of quantum discord close to quantum phase transitions. In contrast to 
the two-spin entanglement, pairwise discord is effectively long-ranged in critical regimes. 
Besides the features of quantum phase transitions, quantum discord is especially feasible 
to explore the factorization phenomenon, giving rise to nontrivial ground classical states 
in quantum systems. The effects of spontaneous symmetry breaking are also discussed 
as well as the identification of quantum critical points through correlation witnesses. 
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1. Introduction 

Quantum mechanics relies on properties of vector spaces, whose the wave function is 
an clement of. Therefore any correlation exploiting the quantum superposition prin- 
ciple (expressing the notion of sum in vector spaces) is indeed of quantum nature. 
This simple evidence has been exploited as a guiding principle to quantify quantum 
correlations beyond the generic notion of 'correlation' in a condensed matter sys- 
tem. In this context, quantum spin systems naturally provide quantum correlations, 
especially at T = 0, where the system is in its ground state. Such correlations are 
diversely important throughout the phase diagram of the system. For example, en- 
tanglement, the most famous quant um c orrelation, turned out strong enough close 
to quantum phase transitions (QPT) Du t it can be very small deep in the phase 
of a highly correlated qua ntum spin system where quantum fluctuations provide 
a factorized ground state (see also Ref. [6] for a review on entanglement in 
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many-body systems). 

A recently devised measure of quantum correlation is the quantum discord 
(QD) 13, providing a 'quantum/class ical s ieve' explicitely based on the superposi- 
tion principle of quantum mechanics I 7 | 8 | 9 [ With QD, it was evidenced that indeed 
there are quantum correlations in separable mixed states. This, in turn, may dis- 
close new scenarios for the phase diagram of many-body quantum systems. The 
first computations of QD for s pin s ystems in the the rmodynamical limit have been 
carried out by Dillenschneider and Sarandy where QD between two spins 
in a mixed state realized by the rest of the system was considered. It turned out 
that the non-analyticities of energy der ivative s at QPTs were reflected not only in 
the two-spin entanglement anc j qd | 10 | 1 1 | j-,^ a j so m |- ne classical correlations, 
establishing them as useful tools to analyze quantum critical phenomena. Furhcr 
inves tigations were then performe d in a number of many-body systems (see, e.g., 
Refs. l i3 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 2 JJ22J) ; where entanglement, QD and classical corre- 
lations were compared. The bulk of evidence we can rely on so far suggests that 
genuine quantum correlations, including entanglement and QD, exhibit the behavior 
expected by scaling theory for the long-range physics close to QPTs. For the short- 
range physics, in contrast, a rather complex situation emerges because entangle- 
ment, QD and standard correlation functions may vanish for different ranges B3 | 24 | 
E]. This is clearly evidenced at factorizing fields occurring in the ordered phase of the 
system (characterized by a nonvanishing spectroscopic gap) . Close to factorization, 
the entanglement range diverges despite the ordinary correlation functions decay 
exponentially^]. 

In this article we review the properties of QD throughout the phase diagram 
of spin systems in a one-dimensional lattice, with nearest-neighbor exchange inter- 
actions. For these systems, different ground states may be achieved through order- 
disorder QPTs occurring via the spontaneous symmetry breaking (SSB) mechanism. 
We consider both the so called thermal states and the symmetry broken ground 
states, enjoying and breaking the symmetries of the Hamiltonian, respectively. As 
long as the system is studied through its energy spectrum, all thermodynamical 
observables stay unaltered irrespective whether the thermal ground state or the 
state with broken symmetry are considered. Nevertheless, the difference comes in 
entanglement and QD (the same argument can be applied to any other measure of 
quantum correlation based on red uced density matrix) . Entanglement was analyzed 
in the SSB scenario in the Refs. | 32 | 33 | 3 4 | 35 | For thc ana i ysis f q Dj SSB effects 
have been first accounted in Refs. I 36 | 3T | -p ne ar tj c i e j s organized as follows. In sec- 
tion [2] we introduce QD. The models we deal with are introduced in the section [3l 
in there we provide a summary of the symmetry breaking mechanism to make the 
article self contained. In the sections |4] and [5] we review the properties of the QD for 



a This, in turn, might be important at low temperature, because indeed thermal fluctuations act 
as infrared cutoff in the system (.the QD clo se t o QPT at low temperature is reviewed in Ref. 1251 
of this focus issue; see also Refs. I 26 l 27 l 28 l 29 l 30 l). 
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thermal and symmetry broken ground states respectively. As i llus tration, in section 
[6] we discuss a witness for multipartite quantum correlations ^H, for both thermal 
and symmetry broken ground states. In section [7] we draw our conclusions. 

2. Quantum discord 

Let us begin by considering a composite system in a bipartite Hilbert space % = 
"Ha ®'Hb- The system is characterized by quantum states described by density 
operators p € B(H), where B(H) is the set of bound, positive-semidefinite operators 
acting on H with trace given by Tr p = 1. Then, denoting by p the density matrix 
of the composite system AB and by p A and p B the density matrices of parts A and 
B, respectively, the total (classical + quantum) correlations between parts A and 
B can be provided by the quantum mutual information 

/(/ : p B ) = S(p A ) - S(p A \p B ), (1) 

where S(p A ) = — Trp A \ogp A is the von Neumann entropy for subsystem A (with 
the symbol log denoting logarithm at base 2) and 

S( P A \p B ) = S(p) S( P B ) (2) 

is the quantum conditional entropy for A given the part B. A remarkable observation 
realized in Ref. [7] is that the conditional entropy can be introduced by a different 
approach, which yields a result in the quantum case that differs from Eq. ^ . Indeed, 
let us consider a measurement performed locally only on part B. This measurement 
can be described by a set of projectors {Bk}. The state of the quantum system, 
conditioned on the measurement of the outcome labelled by k, becomes 

Pk = — {I®B k )p{I®B k ), (3) 
Pk 

where p k = Tr[(J (g) Bk)p{I ® Bk)] denotes the probability of obtaining the outcome 
k and I denotes the identity operator for the subsystem A. The conditional density 
operator given by Eq. (jH]) allows for the following measurement-based definition of 
the quantum conditional entropy: 

S( P \{B k }) = J2PkS(pk). (4) 

k 

Therefore, the quantum mutual information can also be alternatively defined by 

J(p:{B k })=S(p A )-S(p\{B k }). (5) 

Eqs. (TTJ) and ([5]) arc classically equivalent but they are different in the quantum case. 
The difference between them is due to quantum effects on the correlation between 
parts A and B and provides a measure for the quantumness of the correlation. 
Indeed, by generalizing the measure set {B k } to a positive-operator valued measure 
(POVM) and maximizing over such POVMs, we can define the classical correlation 
between parts A and B as^ 

C(p)=m a xJ(p:{B k }), (6) 

{Bk} 
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Then, quantum correlations can be accounted by subtracting Eq. ^ from Eq. ([T]) , 
which yields 

Q(p) = min [l(p A : p B ) - J(p : {B k })} . (7) 

Quantum correlations as measured by Q(p) define QD. As provided by Eq. ([7]), it 
is an asymmetric quantity that measures the quantum correlations between sub- 
systems A and B revealed by measurements on part B. Then, if QD vanishes, the 
system is in a quantum- classical state. In order to define classical-classical states , 
we can symmetrize QD with respect to measurements on each subsystem I 40 | 41 | 
requiring the vanishing of QD for having total classicality. 

Here we apply the above notions to the situation where A and B are individual 
spins of a quantum spin chain. To compute QD, we will adopt the specific case of 
von Neumann local measurements, which arc provided by the orthogonal projectors 
Bk, such that 

B k = Vtl k V^ , fc = 0,l (8) 

where {Ilk} = {\k)(k\} define the computational basis and V € SU(2) is a rotation 
operator. Then, V can be parametrized as 

/ cos f - sin § e"^\ 

where 9 € [0,7r] and cj> € [0, 27 r) are respectively the azimuthal and polar axes of a 
qubit over in the Bloch sphere llQilH. Hence, QD as given by Eq. can be directly 
obtained by minimizing over all angles 8 and <fi. For two spin- 1/2 systems, it is 
shown that the op timal P OVM in Eq. ([7]) is a projective measurement (4-projector- 
elements POVMs) SHH. Remarkably (see Refs. | 44 | 45 | 4 B^ for the part i C ular case of 
Bell-diagonal states (marginal states are maximally mixed), such optimal POVM 
can be achieved through von Neumann measurements, while for the less restricted 
case of ^-symmetric states (X-states), there are examples where optimization over 
POVMs may be inequivalent to von Neumann measurements. In such a situation 
and also for generic nonsymmetric states, which is the case of states with SSB, 
orthogonal measurements arc not more effective than 4-elements POVMs S3 with 
von Neumann measurements in Eq. ([3j providing just upper bounds for QD (as an 
illustration, sec Ref. l48l) . 



3. The models 

Magnetic interactions constitute a simple mechanism to produce classical and quan- 
tum correlations. In this context, we consider the XYZ chain, which is com- 
posed by spin- 1/2 particles localized in a one dimensional lattice, interact- 
ing anisotropically along the three spatial directions with Hciscnbcrg interaction 
strengths J x , J y , J Zl and subjected to a uniform external field h, which is governed 
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by the Hamiltonian 

n =l E [ J * s * s * + J v s i s " + J * s ? s t] - h J2 s * > ( 10 ) 

where Sf (a = x, y, z) are spin-1/2 operators defined on the i-th site of the chain. 
Factorization of the ground state and QPTs characterize the phase diagram of the 
system. The bulk of evidence at disposal so far indicates that the two phenomena 
are interconnected, with the factorization being a precursor of QPTs 



3.1. QPTs and symmetry breaking 

QPTs occur when an external perturbation (like pressure, magnetic field, etc.) 
causes quantum fluctuations that are effective enough to lead to a qualitative change 
of the ground state of the system. In the present article, we consider order-disorder 
QPTs within the symmetry breaking mechanism. Accordingly, the ordered phase 
of the spin system is characterized by a finite magnetization, that is a local order 
parameter breaking the symmetry of the Hamiltonian. In the case of Eq. (|10p , the 
symmetry is a global Zi symmetry generated by rotations of 7r around the z direc- 
tion of all spins: a] -> -<rj = IlVf II, o] -> -a) = ITVJII , where II = (g)^ 1 erf. 
This indicates that the ground states corresponding to opposite {gs\a x ^ v \gs) are 
indeed degenerate and span a ground-state-manifold where any state is a possi- 
ble ground state of the system. Because of such a symmetry (together with the 
fact that H\gs) = e l ^\gs)), the magnetization in the xy vanishes identically in the 
ground state: 

(gs\af y \gs) = (gs\tfa*> y IL\gs) = -<g8\o»> v \gs) . (11) 

The so called 'thermal ground state' is the ground state that is an equal mixture 
of the two degenerate ground states; it can be thought as the limit T — > of the 
thcrmodynamical state. Thermal states enjoy the same symmetry of the Hamilto- 
nian. Superpositions of degenerate ground states may also preserve the symmetry. 
However, they are demonstrat ed t o be not stable to small perturbations because of 
the lack of clustering property S3. Therefore the symmetry is spontaneously broken 
in physical systems (in the thcrmodynamical limit) because any small in-plane local 
magnetic field lifts the degeneracy and there is no observable connecting the two 
ground states (super-selection rule^j). For finite systems, for example, we could 
have a superposition of all the spins in opposed directions like | yt ••• t) i | i~l ■■■ D- 
These two states are orthogonal and preserve the symmetry, but they are not stable 
against small perturbations cither (sec also Rcf . [51]) . Here we will consider both the 
thermal ground states and the state with explicit symmetry breaking of the Zi sym- 
metry. We observe that thermal ground state and the state with broken symmetry 
provide the same thermodynamics of the system (see ^2 for a recent reference); the 
corresponding density matrices, however, differ sensibly. 
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3.2. Factorization of the ground state 



The ground state of Eq.(fT0|) is, in general, a highly entangled state. Nonetheless, 
there may exist points where it is indeed a product state, which is factorized as the 
tensor product of individual spin states For the model Hamiltonian in Eq. (fT0| 
the factorization occurs at 



h f = 2^(J x -J z ){J y -J z ). (13) 

Indeed, a qualitative change of the entanglement, called Entanglement Transition 
(ET) occur s at hf, involving two gapped phases and with a divergent entanglement 
range ESlMI We also mention that analysis of the ground state for finite size systems 
demonstrated that the factorization can be viewed as transition between ground 
states of different parities (see Ref. ^Hl anc j ariC | a related article in this special issue 
. Close to the factorization, the system is characterized by the following pattern 
of correlations functions: 

g a -P(r) = A a -P + B a ^{r){h - h f ) + 0[(h - h f ) 2 } , (14) 

where we remark that A a '" is the same constant V?*. Eq. flTJ| was obtained for the 
quantum XY model (case II below) in Ref. \57\ but it holds for the generic model in 
Eq.([TnD^Sl. 

We will discuss Eq. (fl"U|) in the following cases. 

(I) The non-integrable antiferromagnetic XYX model ( J x = J z = 1), with J y = 
1/4 (this is the case experimentally realized with Cs2 CoCU^^). At zero field h the 
system is critical in the same universality class of the isotropic XY model (see Fig. [T] 
for | J x \ = 1). At finite h the system acquires a finite gap, vanishing at h ~ 3.21 were 
a QPT occurs of the Ising type. The order parameter is (S x ). The factorization is 
displayed at h ~ 3.16. 

(II) For J z — J x = J(l+~f), J y = J(l— 7), J z = 0, the (quantum anisotropic XY) 
Hamiltonian can be diagonalized by first applying th e Jordan- Wigner transforma- 
tion and then performing a Bogoliubov transformation! 6 " ! 61 ! 62 ! The quantum Ising 
model corresponds to 7 = 1 while the (isotropic) XX-model is recovered for 7 = 0. 
In the latter (isotropic) case the model enjoys an additional symmetry resulting in 
the conservation of the total magnetization along the z-axis. The properties of the 
Hamiltonian are governed by the dimensionless coupling constant h/ J (where J is 



b It is remarkable that the factorization occurs for d-dimensional spin system on a bipartite lattice 
and for finite range exchange interaction even in presence of frustration HEl as the field h is varied 
across the factorizing field hf, yielding 

hf = \yJ{J x -J s ){.J y -J z ), (12) 

where J X}V = Z \i-j\ 4V , and J x , y = Z \i-j\ J ^ < Z r indicating the 

range of interaction among the spins. 
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XYX (h=0) 



A' \ 



Fig. 1. Phase diagram of the integrable XYZ model (Eq.JToJ in zero field h = 0). The dashed 
(red) lines correspond to critical phases. Everywhere else the system is gapped. In the points 
A, A' , B, B' the system is tricritical. 

set as the energy scal e and working in units such that h = 1 ) . The phase diagram is 
sketched in Fig.[2p21. In the interval < 7 < 1 the system undergoes a second order 
quantum phase transition at the critical value h c = 1. The order parameter is the 
magnetization in x-direction, (S x ), different from zero for h < 1. In the phase with 
broken symmetry the ground state has a two-fold degeneracy reflecting a global 
phase flip symmetry of the system. The magnetization along the z-direction, (S z ), 
is different from zero for any value of h, but has a singular behavior in its first 
derivative at the transition point. In the whole interval < 7 < 1 the transition 
belongs to the Ising universality class. For 7 = the QPT is of the Berezinskii- 
Kosterlitz-Thouless type. The ET transition, occurring at the the factorization field 
is also associated with a change in the way the correlations decay with dist anc e: for 
h < hj, the decay is oscillatory while, for h > hf, the decay is monotonic ^21. The 
quantum XY model can be experimentally realized with the magnetic compound 
CoNb 2 6 E1. 



(Ill) The antifcrromagnctic anisotropic XX Z chain can be obtained from Eq. (|10| 
by setting , J x = J y = J, and J z = A. At zero field it presents a critical xy phase 
with quasi-long range order (quasi-lro) for \J Z \ < 1 ; this is separated by two classi- 
cal phases with QPTs at J z = ±1. For h ^ the xy phase is a strip in the phase 
diagram, eventually turning into polarized phases for sufficiently strong magnetic 
field. Here we remark that the factorization phenomenon degenerates in the satu- 
ration occurring as a first order transition. Despite the similarities the factorization 
and saturation are quite distinct phenomena. In the case in which hj J = and 
for any value of A the model is referred to as the XX Z model. The two isotropic 
points A = 1 and A = — 1 describe the antiferromagnetic and ferromagnetic chains 
respectively (see the phase diagram in Fig.[3p^|). The one-dimensional XX Z model 
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Fig. 2. The zero temperature phase diagram of the one dimensional anisotropic XY model in 
transverse field. Along the quantum critical line the model identifies the Ising universality class 
with indices z — v = 1; in the hatched area (with 7 > 0) the system display long range order in 
the x — y spin components. The critical XY regime coicide with that one of the XXZ model for 



A = 0. The factorization of the ground state occurs along the circle h = \/~\~- 



■T 



can be solved exactly by Bethe Ansatz technique (see e.g. ^2) anc j the correlation 
functions can be expressed in terms of certain determinants (see for a review). 
Correlation functions, especially for intermediate distances, are in gene ral diffi cult 
to evaluate, although important steps in this direction have been made EElfiZI 




Fig. 3. Zero temperature phase diagram of the spin 1/2 XXZ model in one dimension. The 
XY phase is characterized by power law decay of the xy correlations. The Noel and the XY 
phases are separated by a line of second order phase transitions; at A = 1 the transition is 
of Berezinskii-Kostcrlitz-Thouless. The ferromagnetic and XY phases are separated by first order 
phase transitions due to simple level crossings; the onset to the ferromagnetic phase occurs through 
the saturation phenomenon. 



The zero temperature phase diagram of the XXZ model at zero magnetic field 
shows a gapless p hase i n the interval — 1 < A < 1 with power law decaying cor- 
relation functions I 68 | 69 [ Outside this interval the excitations are gapped. The two 
phases are separated by a Bcrczinskii-Kostcrlitz-Thoulcss phase transition at A = 1 
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while at A = — 1 the transition is of the first order. In the presence of the external 
magnetic field a finite energy gap appears in the spectrum. The universality class 
of the transition is not affected, as a result of the conservation of the total spin-z 
component ^31. 



3.3. The reduced density matrix 

In order to compute the QD Q r between any two spins A and B at distance r 
along the chain, the key ingredients are the single-site density matrices p ,p B and 
the two-site density matrix of the composite subsystem p AB [see Eq.Q]. Due to 
translational invariance along the chain, single-site density matrices are the same 
for any spin and they are given by 

~A ~B 1 ( 1 + 9z 9x \ n K , 

where g a = (<r Q ) are the local expectation values of the magnetization along the 
three different axes. 

The two-site reduced density matrix for a Hamiltonian model Eq. (|10j) reads 



Pr = i 



/A a a F\ 

a B C b 

a C B b 

\F b b D/ 



(16) 



in the basis {|00), |01), |10), |11)}, where |0) and |1) are eigenstates of a z (because of 
translational invariance, this density matrix depends only on the distance r between 
p AB ). The various entries in Eq. (fl~6|) are related to the two- 



9zz, 
9zz, 



the two spins: p r = 
point correlators gap{r) 
the following: 

A = 1 + 9 
D = 1 - 9 
B = 1 - g zz 
C 9xx ~t~ 9yyi 
E 9xx 9yy- 
The entries above express the parity coefficients, while 

a = 9x + 9xz j 
b = g x - 9xz ■ 



and to the local magnetizations, according to 



(17) 



(18) 



As long as the state displays the ^-symmetry of the Hamiltonian (i.e., the thermal 
ground state) the matrix elements a and b vanish. In this case, the reduced density 
matrix correspond to an X state. On the oth er han d, the requirement a,b =/= 
characterizes the state with symmetry breaking l ^4 | 35 | 
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4. Quantum discord in thermal ground states 

In this section we summarize the properties of QD in the thermal ground states 
of spin models with Z 2 -symmetry. In this case the density matrices are given by 
requiring a = b = into Eq. ()16[) . 

4.1. XY model 

We begin with the transverse field XY chain (case (I) of Sec. |3]). Here, we will be 
interested in the correlations between arbitrary distant spin pairs (not only nearest- 
nighbor pairs). Then, the reduced state po r for spin pairs at a distance r = \i — j\ 
reads 

P0r = |{I° n + («T Z )K+<)+ E <^>0<}, (19) 

where I 0ra is the identity operator acting on the joint state space of the spins and 
n. The magnetization density (a z ) as well as the two-point correla tion functions 
((To<7*) can be directly obtained from the exact solution of the model ^1 Then, the 
total information shared by the spins in the state (|19|) is given by 

I(p 0r ) = S(p ) + S(p r )-S( Pr ), (20) 



with 



and 



where 



q/ \ qf \ ^ 1 + (-1) V) , l + (-l)V) 

S(po) = S(p r ) = -2_^ 7, fog 2 S (21) 



i=0 



S(por) = lo S2 & + Vi lo S2 Vi), (22) 



6 = 1 + ( f^ + ^^/(( ff0 X)-W» 2 +4( ff f, 

% = \ [1 - + (-!)*(«<) + W»] • (23) 

Following | H | 24 | 14 l^ we obtain that the minimum in Eq. ([7j is attained, for all 
values of h, 7, and n, by the following set of projectors: {|+)(+|, |— )(— |}, with 
|±) = (|t) ± \l))/\/2, where {|f) , \l)} are the eigenstates of o z . Thus one obtains 

J {fir) = Hunifil) + H bin (p 2 ), (24) 

where H b i n {x) is the binary entropy 

H bln (x) = -x\og 2 x- (1 - x)log 2 (l - x) (25) 

and 

Pi = \{l + (<r z )), P2 = ^l + V / K T ^) 2 + (^) 2 )- (26) 
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This provides therefore an analytical expression for evaluating QD. We can then 
use this expression to investigate QD at quantum criticality. By fixing 7 = 1 (Ising 
model) and considering pairwisc QD between first neighbors, we can identify the 
ferromagnetic-paramagnetic QPT by looking at cither classical or quantum corre- 
lations, as shown in Figs. [4] and [5] 




Fig. 4. (Color online) First derivative ol the classical correlation lor nearest-neighbor spins with 
respect to h in the transverse field Ising chain lor different lattice sizes L. The derivative ol C 
has a pronounced minimum at h m i n , which tends to the critical point h = 1 as L — > oo. Inset: 
dC/dh taken at h m i n exhibits a logarithmic divergence fitted by (dC/dh)\ h . = —0.29161 — 
0.22471 log L. [Original plot from Ref. 11] 




h 



Fig. 5. (Color online) First derivative ol the quantum correlation lor nearest-neighbor spins with 
respect to h in the transverse field Ising chain lor different lattice sizes L. Inset: dQ/dh presents 
an inflection point that tends to the QCP h = 1 as L — > oo. [Original plot from Rel. 11] 



As easily observed, derivatives of both C(p) and Q(p) exhibit a non-analytical 
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behavior at the quantum critical point as we approach to the limit of an infinite 
system. 
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Fig. 6. Decay of quantum discord with the distance between the spins sites for some values of h 
and 7. The points are the computed values of QD and the lines are the exponential fits (see the 
text for details). [Original plot from Ref. 24] 
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Remarkably, besides the non-analytical behav ior of pairwise QD derivatives, we 
can show that QD exhibits a long-range decay This is in contrast with the 
two-spin entanglement behavior, which is typically very short-ranged. In order to 
illustrate the long-range behavior of QD, we plot in Fig. [6] pairwise QD as a function 
of distance n. The curves in Fig. [6] are the exponential fits of QD. We observe that, 
for both examples of anisotropics considered, 7 = 0.1 and 7 = 0.5, the decay of QD 
with distance can be well fitted by an exponential function a + 6cxp(— cr), where 
a, b, and c are constants. Nevertheless, we notice that, while for h > 1 QD vanish 
exponentially, in the cases where h < 1, we obtain a constant long-distance value 
for QD that depends only on 7 and h. It is also remarkable to observe that the 
factorization phenomenon can be traced already for the thermal ground state: it is 
the unique value of the field where the same quantum correlations are present at any 
length scale (left inset of Fig. EJ. This is due by the peculiar pattern of correlation 
functions close to the factorizing field in Eq. (fl4|) . 



4.2. XXZ model 

As a second example, let us consider now XXZ spin- 1/2 chain with h = (the 
case (III) of Sccj3]). In this case the Hamiltonian exhibits U(l) invariance, which 
ensures that the element F of the reduced density matrix given by Eq. (fT6|) vanishes. 
Moreover, the ground state has magnetization density g k z = (crjf) = (Vfc), which 
implies that A = D = (1/4) (1 + g zz ). For convenience, we consider correlations 
between nearest-neighbor spin pairs, i.e. p r = p\ = p, and then expand p in terms 
of Pauli operators, which reads 



i=l 



(27) 



with 



ci = c 2 = - (g xx + g yy ) , c 3 = g z 



(28) 



By taking into account those notations, classical correlations in Eq. © can be 
analytically worked out yielding 



C{ P ) = log (1 - c) + log (1 + c) , 

with c = max(|ci|, | C2 1 , | C3 1 ) . For the mutual information I(p) we obtain 



(29) 



3 

7(p)-2 + ^A 4 logAi, 

i=0 



(30) 
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where 

A = - (1 - Cl - C2 - c 3 ) , 
Ai = A 2 = i (1 + c 3 ) , 

A 3 = \ (1 + d + c 2 - c 3 ) . (31) 

In order to compute C(p) and Q(/o) we write ci, C2, and c 3 in t erms o f the ground 
state energy density. By using the Hellmann-Feynman theorem l ^l l ^ l for the XXZ 
Hamiltonian Hxxz, we obtain 

1 de xxz 
c i = c 2 = 2 (G X x + G yy ) = A e^j-j , 

c 3 = G« = -2^, (32) 



<9A 

where e xxz is the ground state energy density 

£XXZ = ■ ^ == _ 2 ( XX + Cyy + AG ZZ ) , (33) 

with |-0o) denoting the ground state of Hxxz- Eqs. (p?2"j) and ([3"5| hold for a chain 
with an arbitrary number of sites, allowing the discussion of correlations either for 
finite or infinite chains. Indeed, ground state energ y a s well as its derivatives can 
be exactly determined by Bethe Ansatz technique which allows us to obtain 
the correlation functions c\, c 2 , and c 3 . In Fig. we plot classical and quantum 
correlations between nearest-neighbor pairs for an infinite XXZ spin chain. 




Note that, in the classical Ising limit A — > oo, we have a fully polarized ferro- 
magnet. The ground state is then a doublet given by the vectors | tt " " " t) an d 
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| XX ■ ■ ■ l) , yielding the mixed state 

p = \\ tt •■■ tXtt •■■ 1 1 + \\XX ■ ■■ X)(U ■■■ I I- (34) 

Indeed, this is simply a classical probability mixing, with C(p) = I(p) = 1 and 
Q{p) = 0. The same applies for the antifcrromagnctic Ising limit A — > — oo, where a 
doubly degenerate ground state arises. Moreover, observe that the classical (quan- 
tum) correlation is a minimum (maximum) at the infinite-order QCP A = — 1. On 
the other hand, both correlations are discontinuous at the first-order QCP A = 1. 
The behavior of QD at the first-order QPT is in agreement with the expected be- 
havior for entanglement, as shown by Wu et al. in Ref. 1121 In fact, the approach 
used by Wu et al. for entanglement also applies for any function of the reduced 
density matrix of two-spins. It is only based on the fact that the non-analyticitics 
in the derivatives of the ground state energy typically come from the elements of 
the reduced density matrix and, as such, may be inherited by entanglement or other 
functions of the density matrix elements. There is a caveat in the fact that entangle- 
ment and QD definitions involve a maximization or minimization processes, which 
may create accidental non-analyticities or even hid e ones. In fac t, th e behavior of 
the entanglement as measured by concurrence ^1 or negativity^] at A = 1, is 
not the one expected for first-order QPTs, but rather for second order QPTs, since 
it is not discontinuous at A = 1 but its derivative is. Such a misidentification of 
the order of the QPT is due to the maximization process involved^! Remarkably, 
although involving an optimization, QD does indicate the correct order of the QPT, 
exhibiting in this case a behavior that is superior to entanglement. For A = — 1, the 
result for QD is in agreement with the entanglement behavior (maximum or mini- 
mum at the infinity-order QPT), even though there is no proof of the generality of 
this behavior for arbitrary infinite-order QPTs. 

5. Quantum discord in ground states with symmetry breaking 

In this section, we discuss the main features of QD for p r given by Eq. (|16[) with 
non vanishing a, b given by (fT8|) . This correspond to the ground state of the Hamil- 
tonian (|10j) (with nearest neighbor couplings), with broken symmetry. In this case 
the reduced density matrix is accessed via DMRG for finite systems with open 
boundaries by adding a small symmetry-breaking longitudinal field (of typi- 
cal magnitude h x = 10~ 6 ) to the Hamiltonian fl In Fig. it is displayed QD for 
the quantum XY model. Similarly, QD for XYX model in the ground state with 
symmetry breaking is obtained in Fig.®. 

QD is substantially affected by SSB We observe that quantum correlations 
are typically much smaller deep in the ordered ferromagnetic phase h < h c , rather 



c For the quantum Ising model the reduced density operator with o,l) / could be accessed 
analytically, in principle. Nevertheless, we resort to DMRG because the g X zix) is given in terms 
of integral contour formulas E51. 

d In contrast, bipartite entanglement is not affected by SSB around the quantum critical point. SSB 
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h 

Fig. 8. Quantum discord Q r (h) between two spins at distance r in the XY model at 7 = 0.7 
(main plot and left inset) and 7 = 1 (right inset), as a function of the field h. Continuous lines 
are for the thermal ground state, while symbols denote the symmetry-broken state obtained with 
DMRG in a chain of L = 400 spins; simulations were performed by keeping m = 500 states and 
evaluating correlators at the center of the open-bounded chain. For 7 = 0.7 and at hf ~ 0.714, 
in the symmetric state all the curves for different values of r intersect, while after breaking the 
symmetry Q r is rigorously zero. At the critical point Q r is non analytic, thus signaling the QPT. 
In the paramagnetic phase, there is no symmetry breaking to affect Q r . [Original plot from Ref. 36] 




h 

Fig. 9. Quantum discord Q r (h) between two spins at distance r in the XYX model. 

than in the paramagnetic one h > h c . Nonetheless, as we shall see, they play a fun- 
damental role to drive the order-disorder transition at the QPT, where Q r exhibits 



slightly influences bipartite entanglement only for h < hj. However, multipartite entanglement is 
highly affected ESI. 
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Fig. 10. (a) Finite-size scaling of dhQi for the symmetry-broken state in proximity of the critical 
point h c . Displayed data are for 7 = 0.7. The first derivative of the QD is a function of L~ v (h—hm) 
only, and satisfies the scaling ansatz duQl ~ x F[L~ v (h — h m )] , where /i m is the renormalizcd 
critical point at finite size L and ui = 0.472. We found a universal behavior h c — h m ~ 1.28±0.03 
with respect to 7. Inset: raw data of dhQl as a function of the transverse field, (b) Scaling 
of Qi close to the factorizing field, for 7 = 0.7: we found an exponential convergence to the 
thermodynamic limit, with a universal behavior according to e~ aL (h — ), a R3 1 [hj denotes 

the effective factorizing field at size L, while <5(Qi) = Q\ L ^ — Q' i_> °°^]. Due to the extremely fast 
convergence to the asymptotic value, already at L ~ 20 differences with the thermodynamic limit 
are comparable with DMRG accuracy. Inset: raw data of Qi as a function of h. The cyan line is 
for L = 30 so that, up to numerical precision, the system behaves at the thermodynamic limit. 
[Original plot from Ref. 36] 



a maximum, as well as the correlation transition at hf, where Q r is rigorousl y zero. 
The QPT is marked by a divergent derivative of the QD (see also l 1 ^ ! 11 ! 14 !). Such 
divergence is present at every 7, for the symmetry broken state; on the other hand, 
for the thermal ground state, it is not present at 7 = 1. A thorough finite-size scaling 
analysis is shown in Fig. (fTUT b)) proving that z = v = 1. For the thermal ground 
state (in the thermodynamic limit), we found that dhQ r diverge logarithmically as 
dhQr ~ In \h — h c \, within the Ising universality class. 
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At the factorizing field hf, all the correlation measures are zero in the state 
with broken symmetry (see symbols in Fig. [HJ ; in particular, we numerically found 
a dependence 

Q r ~(h-h f ) 2 x(^L) r (35) 

close to it. Such behavior is consistent with the expression of correlation functions 
close to the factorizing line Eq. (fT4|) . and here appears incorporating the effect arising 
from the non vanishing spontaneous magnetization. 

It is found that Q r exponentially to asymptotic value Q^^ 00 ^ (sec Fig. rTUT b^. 
j n | 54 | 53 |^ jj. wag s j 10wn h j marks the transition between two different patterns of 
entanglement. The factorization is thus a new kind of zero-temperature transition of 
collective nature, not accompanied by a change of symmetry. We emphasize, though, 
that this transition docs not correspond to any non analyticity in the ground state as 
a function of h. 0. At the level of spectral properties of the system, we interpret this 
result as an effect o f ce rtain competition between states belonging to different parity 
sectors for finite as these states intersect, the ground-state energy density is 

diverging for all finite L (such divergence, though, vanishes in the thermodynamic 
limit). 



6. Witnessing quantum correlated states 

Computation of QD involves an extremization procedure. However, if we just want 
to find out whether a state exhibits quantum correlations, such a procedure can 
be avoided by introducing witness operators, namely, observables able to detect 
the presence of nonclassical states. A witness W is defined here as a Hermitian 
operator whose norm ||W|| is such that, for all nonclassical states, || W\\ > 0, with a 
convenient norm measure adopted. Therefore, ||W|| > is a sufficient condition for 
nonclassicality (or, equivalcntly, a necessary condition for classicality) . Let us begin 
by defining a classical state. We will concetrate in von Neumann meansurements, 
but our approach can be generalized for POVMs. 

Definition 1. If there exists any measurement {H,} such that &(p) = p then p 
describes a classical state under von Neumann local measurements. 

Therefore, it is always possible to find out a local measurement basis such that a 
classical state p is kept undisturbed. In this case, we will denote p £ C N , where C N 
is the set of iV-partite classical states. Observe that such a definition of classicality 
coincides with the vanishing of global QD EH. Note also that $(p) by itself is a 
classical state for any p, namely, $($(p)) = $(/o). Hence, $(p) can be interpreted 
as a decohered version of p induced by measurement. A witness for nonclassical 
states can be directly obtained from the observation that the elements of the set 



e Accordingly, the fidelity ^(h) = \ (gs(h)\gs(h + Sh)\ (which can detect both symmetry breaking 
and non-symmetry breaking QPTs), is a smooth function at hf. 
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{Uj} arc cigcnprojectors of p. This has been shown in Rcf. [5H being stated by the 
theorem below (see also Ref. [70j) . 

Theorem 2. p e C N ^ [p,Uj] = (Vj), with Uj = U\ <8 - and j 

denoting the index string (ii • ■ ■ in) ■ 

Proof: If p € C N then $(p) = p. Then, similarly as in Ref. [701 a direct evaluation 
of Ej \P, n,l [p, Ilj? yields £\ [p, n,] [p, IL/ = $(p 2 ) - p 2 . However, $(p) = p also 
implies that $(p 2 ) = p 2 . Therefore, 

5>,H,] [p.II^O. (36) 

i 

Hence, p e C w [p, IT,-] = 0. On the other hand, if [p, EL,-] = then {n^} provides 
a basis of eigenprojectors of p. Then, from the spectral decomposition, we obtain 

3 h,— i*iv 

which immediately implies that $(p) = p. Hence [p, Hj] = =>■ p € C w . □ 

We can now propose a necessary condition to be obeyed for arbitrary multipar- 
tite classical states. 

Theorem 3. Let p be a classical state and pA t the reduced density operator for the 
subsystem Ai. Then [p, pa x <8> • ■ ■ <8 Pa n ] = 0. 

Proof: From theorem [21 if p € C N then the spectral decomposition of p yields 
Eq. ([37]). Therefore, 

PA s =J>« n X- (38) 

Hence, by direct evaluation, we obtain that [p, pa x ® ■ ■ • ® PA N ] =0. □ 

Observe that, given a composite multipartite state p, it is rather simple to eval- 
uate the commutator [p, pa 1 <£> • • • ® Pa n ], with no extremization procedure as usu- 
ally required by QD computation. From the necessary condition for classical states 
above, we can define a witness for nonclassicality by the norm of the operator 
[p, pA! ■ ■ ■ PA N ]- Indeed, if p € C N \\W\\ = 0, where 

W= [p, p Al ® ■ ■ ■ <Z> pa n ] ■ (39) 

Therefore, ||W|| > is sufficient for nonclassicality. 

For concreteness, we will take ||W|| as defined by the trace norm, namely, 

|| W|| = TrVWWt . (40) 

In order to apply the witness ||W|| in the XY model, we will first consider the 
case of the thermal (^-symmetric) ground state. For this case, a plot of ||W|| as a 
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function of h is provided in Fig. [IT] for 7 = 0.6. Observe that this plot shows that 
the only possible classical state appears for h = 00. Indeed, this point corresponds 
to an infinite transverse magnetic field applied, which leads the system to a product 
state, with all spins individually pointing in the z direction. 

However, besides large h, the ground state with SSB exhibits a further nontrivial 
factorization (product state) point at -f 2 + h 2 = l.This point identifies a change in 
the behavior of the cor relation functions decay, which pass from monotonically to 
an oscillatory decay ^1 Besides, it has been realized that, at this poin t, the product 
ground state is two- fold degenerated even for finite systems I 79 | 80 | 36 | 



0.25 
0.2 

0.15 

0.1 

0.05 



"0 0.5 1 1.5 2 2.5 

l/h 

Fig. 11. (Color online) Witness ||W|| and exact value of QD for a pair of spins in both symmetric 
and broken ground states of the XY spin-1/2 infinite chain as a function of h for 7 = 0.6. [Original 
plot from Ref. 38] 

On the other hand, if SSB is taken into account, which must be the case in 
the thermodynamic limit, this picture dramatically changes. Indeed, in the broken 
case, the two-spin reduced density matrix is more complicated (see Eqs. (fTB")) . ([18])). 
Tight lower and upper bounds for a{h, 7) and b(h, 7) are obtained in Ref. 1351 In our 
plots, either of such bounds essentially yields the same curves, which led us to keep 
the results produced with the lower bound. Concerning the results for ||W||, we can 
compare the plots for the thermal ground state and the ground state with SSB also 
in Fig. [TT] Observe that, in the case with SSB, the only point for which ||W|| = 
is l/h = 1.25. Indeed, this corresponds to a classical state, as can be confirmed by 
the exact computation of QD. Since our witness identifies classical- classical states, 
we adopted the symmetric version of QD | 40 | 41 | R em arkably, this classical state is 
associated with one of two doubly-degenerated factorized (fully product) ground 
state, which appears as a consequence of the Z2 SSB. This unique classical point 
for the XY model is then clearly revealed by the witness evaluation. Another aspect 
of the witness observed from Fig. [TT] is that it displays nonanalyticity in its first 
derivative at the QPT h = 1. This phenomenon occurs both for the symmetric and 
broken ground states, which promotes ||W|| to a simple tool also to detect QPTs. 



* — a Witness with Z 3 symmetry 
■ " Witness with SSB 

Exact discord with SSB 
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7. Conclusion 

We reviewed on the two-spin QD in the ground state of spin chains of the XYZ type 
in a transverse field. Close to QPTs, the anomalies of QD reflect the universality 
class of the Hamiltonian of the system. Therefore we can conclude that the scal- 
ing of QD and entanglement is identical with the known scaling of correlations of 
standard many body theory. In contrast, the range of quantum correlations at the 
critical point is non universal: entanglement is typically short-ranged (with scaling 
laws depending on the microscopic details of the system), while QD decays alge- 
braically (sim ilarly as the standard correlation functions) following the universality 
1941 

paradigm L ^- ±1 . In particular, we note that QD is generically more robust than entan- 
glement. This ultimately arises because QD catches quantum correlations in separa- 
ble (mixed) state, where entanglement is vanishing. Therefore, at zero temperature, 
the QD can be sensible at distances were entanglement is vanishing. Along simi- 
lar lines of reasoning, QD can survive at temperatures where the states arc mixed 
enough to kill the entanglement (see Ref. [5S]m thi s speci al issue). The investigation 
of multipartite measures of quantum correlations III as well as their impact in 
the characterization of QPTs is certainly a future challenge. Entanglement is be- 
lieved entering in its multipartite form close to the QPT; the two-spin entanglement 
is then a small amou nt of t he entanglement at disposal, ultimately because of the 
monogamy property I 32 | 83 [ Q n other hand, it has alrea dy been shown that QD 
does not obey the usual monogamy relationship | 84 | 85 | 86 [ However, a generalized 
monogamy constraint can be established^! The implications of such monogamy for 
the behavior of multipartite QD at QPTs is object of current research. For studies 
of multipartite QD in spin models, see Refs. [41] [82j Moreover, the application of 
genuine multipartite measures of QD (e.g., as defined by Refs. [87, 88, 89) is also a 
further promising topic. 

At the factorizing field, the correlations do not depend on the spatial distance. 
Therefore all the QDs for different ranges crosses at h = hf (see the inset of FigEJ). 
The actual value of the discords is finite for thermal ground states and it is van- 
ishing for ground states with symmetry breaking. This results since the factorized 
thermal states are mixed; in the case of symmetry breaking the factorized two-spin 
state is pure, and therefore QD coincides with entanglement (they are both vanish- 
ing at hf). Here we remark that QD is a smooth function close to the factorizing 
field (in contrast, entanglement displays certain discontinuity), and therefore it can 
be analyzed easly (see FigllObl ^^^^ l In fact the finite size scaling of the factor- 
ization phenomenon was achieved through the QD and not through the analysis 
of entanglement. We finally comment that the role of multipartite entanglement 
was demonstrated to be negligible close to the factorizing field ^31. Therefore the 
two-spin entanglement developed so far can provide a complete enough scenario of 
entanglement close to hf. 

To conclude, we mention that QD is also interesting from the experimental 
perspective. Pairwise non-classical correlations in highly mixed states were exper- 
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imentally measured in NMR systems In that experimental protocol non linear 
witness operators (see section [5]) were employed to simplify the optimization en- 
coded in the discord. Interestingly enough the QD can be easly evaluated for spin 

cluster. This may be important to measure the non-classical correlation via neutron 
|qo| 

scattering- 2 ^ 1 . 
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